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ABSTRACT 

Recent ultraviolet and X-ray observations of intergalactic OVI and OVII absorption systems 
along lines of sight to bright quasars have opened a new window onto the "warm-hot in- 
tergalactic medium" (WHIM). These systems appear to provide a significant reservoir for 
baryons in the local universe, and comparison to cosmological simulations suggests that their 
abundance roughly matches theoretical predictions. Here we use analytic arguments to elu- 
cidate the physical properties of the absorbers and their role in structure formation. We first 
show that if the absorbers result from structure-formation shocks, the observed column densi- 
ties naturally follow from postshock-cooling models, if we include fast-cooling shocks as well 
as those that cannot cool within a Hubble time. In this case, the known OVI absorbers should 
show stronger OVII absorption than expected from collisional-ionization equilibrium (and 
much more than expected for photoionized systems). We then argue that higher-temperature 
shocks will be spatially associated with more massive virialized objects even well outside the 
virial radius. Thus the different oxygen ions will trace different structures; OVII absorbers are 
the most common because that ion dominates over a wide temperature range (corresponding 
to a large range in halo mass). If each dark-matter halo is surrounded by a network of shocks 
with total cross section a few times the size of the virialized systems, then we can reproduce 
the observed number densities of absorbers with plausible parameters. A simple comparison 
with simulations shows that these assumptions are reasonable, although the actual distribution 
of shocked gas is too complex for analytic models to describe fully. Our models suggest that 
these absorbers cannot be explained as a single-temperature phase. 
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1 INTRODUCTION 

The cosmic-web paradigm has had a great deal of success in ex- 
plaining the distribution of matter in the intergalactic medium 
(IGM) at moderate redshifts z ~ 3, especially in the Lya forest 
(e.g.. iRauchll 1 998^ . In this picture, dark matter collapses through 
gravitational instability into sheets and filaments before accreting 
onto bound structures such as galaxies, groups, and clusters (which 
appear at the intersections of filaments). The baryons also accrete 
onto these structures, but (in addition to gravitational forces) they 
are subject to fluid effects such as shock heating. At z ~ 3, the 
infall velocities are modest relative to the ambient sound speed 
and shocks are relatively unimportant except in collapsed objects. 
However, at the present day, shocks h ave heate d mo st of the fil- 
amentary gas to T>10 5 K JCen & Ostriked fl999): this gas is 
now known as the "warm-hot IGM" (WHIM). Simulations pre- 
dict that this phase contains a su bstantial fraction of the baryons 
Jcen & Ostrikei 1999; Dave et al. 2001). The characteristic tem- 
perature T ~ 10 K corresponds to the postshock temperature that 
has collapsed over the nonlinear mass scale in one Hubble time. 
Unfortunately, most of the WHIM lies at a modest density 
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(~ 10-100 times the cosmic mean) and in a temperature regime 
with few observational signatures; the hydrogen column densi- 
ties are small because of collisional ionization and the gas is not 
hot enough to emit substantial bremsstrahlung radiation. As such, 
the WHIM has thus far only been detected through absorption 
of highly-ionized metals (specifically OVI, OVII, and possibly 
OVIII). Space-based ultraviolet spectrographs (including STIS on 
the Hubble Space Telescope and the Far Ultraviolet Spectroscopic 
Explorer) have allowed relatively straightforward detections of the 
OVI AA 1032, 1038 double t along lines of sight to nearby bright 
uasars jTrip p et al] l2000t jTCrjr^^^avageJj2p0p^JOe^erle^^l] 



200(1 ISembach et aimOOlj: Isavaee et alj 



Rich ter et alj|2004 fP rochaska et al. 2004 



2002; Shulletal. 2003 



Tumlinson et al. 2005 



Danforth & Shull 2005). The inferred number density of absorp- 
tion systems is large (dn/dz ~ 15), suggesting that these systems 
constitute a significant baryon reservoir. However, in collisional- 
ionization equilibrium (CIE), OVI only exists near the lower end 
of the WHIM temperature range. According to the simulations, 
most of the gas can only be probed through OVII or OVIII ab- 
sorption, both of which require X-ray studies. The current gener- 
ation of instruments does not have the sensitivity to measure the 
weak absorption expected from the IGM, and there are relatively 
few firm detections. A number of groups have found OVII and 
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OVIII absorptio n at z = that may be associa ted with t he Lo- 
cal G roup iNicastro et alj 2002t iRasmussen et all2003t iFang et all 
120031) . iNicastro et alJ 120051) recently discovered two OVII ab- 
sorbers along the line of sight to the blazar Mkn 421 during an out- 
burst. Other detections have relatively low significance or lack con- 
firmation from different instruments (e.g., IFang etalj2 002). Never- 
theless, the Mkn 421 line of sight suggests an even higher number 
density than OVI (dn/dz ~ 60), albeit with large statistical uncer- 
tainties. 

Comparison to lines of sight extracted from simulations sug- 
gests that these number densities are comparable to theor etical 
expec t ations iHel lsten et aT]|l998j: |Cen et al Jl200 it IFang & Brvanl 
1200 it IFang et alJ 120021: IChen et alJl2003l) . Detailed quantitative 
tests are difficult because of the large uncertainties on both the 
observational and theoretical sides (with the latter depending pri- 
marily on the unknown distribution of metals in the IGM), but 
the overall consistency is a reassuring confirmation of the cosmic- 
web paradigm. However, the existing studies have not addressed 
many questions about their physical origin or characteristic envi- 
ronment. As a result, we lack a clear picture of their specific role 
in the structure-formation process. For example, do they probe out- 
lying filamentary gas, gas more closely associated with galaxies, 
or something entirely different? Do the OVI and OVII absorbers 
probe distinct environments, and do the two ions have different 
physical origins? What sets the observed column densities? 

In this paper, we provide an analytic framework to answer 
these questions. This is clearly a difficult proposition: the WHIM 
is nonlinear and inherently asymmetric, so it is difficult to model 
analytically. Nevertheless, models describin g the gross pro perties 
of the WHIM do exist iNath & Silkll200ll: IValageas et alj|2002t 
IFurlanetto & Loebll2004l) . and it is interesting to consider what 
these simple arguments reveal about the absorbers. Moreover, if 
shocks are ultimately responsible for heating gas to the WHIM tem- 
peratures, simple models of postshoc k cooling can h elp to explain 
the properties of the absorbers iHeckman et all2002h . We will con- 
nect these two approaches in order to identify the key properties of 
the absorbers and how they may relate to collapsed objects. 

The plan of this paper is as follows. In Sj2| we show how post- 
shock cooling predicts characteristic column densities for the ab- 
sorbers. We then construct a simple analytic model to predict the 
abundance of IGM oxygen absorbers in ij3| We highlight some of 
the strengths and weaknesses of this model by comparing to cos- 
mological simulations in Sj4] In f|5] we show that comparing the 
OVI and OVII columns can robustly distinguish shock-heated and 
photoionized OVI absorbers. Finally, we conclude in ij6| 

We will primarily be concerned with OVI, OVII, and OVIII 
in this paper. For reference, we list details of their strongest tran- 
sitions here. OVI has a UV doublet at A = 1032, 1038 A; the 
corresponding oscillator strengths are f^ 1 = 0.133, 0.0660. The 
most prominent OVII line is a Hea transition with A = 21.602 A 
and /2y n = 0.6945, while that for OVIII is a Lya transition with 
A = 18.97 A and f™ u = 0.416. To avoid any assumptions about 
the velocity widths of each line, we will present all of our results in 
terms of column density Mon. Observations are instead often pre- 
sented in terms of equivalent width W (which is the fundamental 
observable). In the unsaturated limit, the two quantities are related 
(in cg s units) via W = 8.85 x 10~ 13 A?ion/o S c A 2 (e.g.. ISpitzeJ 
1978). (For OVI, W is usually, though not universally, quoted for 
the stronger of the two doublet transitions.) 

Throughout our discussion we assume a cosmology with 
Q m = 0.3, Q.A = 0.7, D. b = 0.046, H = WOh km s" 1 Mpc" 1 



(with h — 0.7), n = 1, and ag = 0.9, consistent with the most 
recent measurements fcpergel et all2003h . 



2 THE CHARACTERISTIC COLUMN DENSITIES 

We will begin with the simple question of whether shock heating 
can account for the observed column densities of OVI absorbers. 
A number of recent observations have found a large number of 
such systems in the local universe, with typical column densities 
in the range JVpyi « 2-20 x 10 13 cm" 2 (e.g.jTripp et alfcOOfl 
ISembachetai]|200lL iRichter et al]|2004l : IProchaska et alj|2004h . 
W hile the lower limit may b e an observational selection effect (but 
see Danfort h & Shulll2005l) . the rapid decline in abundance as the 
column density approaches 10 14 cm -2 suggests that this is a char- 
acte ristic maximum strengt h. 

iHeckman et alJ J2002) provide a simple argument for why we 
should expect shocks to produce columns in this range. Consider a 
parcel of gas that is shocked to a temperature T s . It will flow away 
from the shock with a velocity v g (fixed by the jump conditions). 
After a time t coo \ = 3kT s /n s A it will have radiated its thermal 
energy and no longer be highly ionized; here n s is the postshock 
gas density and A its normalized cooling rate in ergs cm 3 s _1 . If 
we view the system along a plane perpendicular to the shock, and 
(/ the gas cools fully, we therefore expect to see highly-ionized gas 
along a path length Z coo i ~ v g t zoo i with a corresponding column 
density N ~ n s l coo i, Thus 

Novi = /ovi Z [O/H]© n s v g i coo i (1) 

where [O/H] = 10" 3 07 is the oxygen abundance by num- 
ber in gas with solar m etallicity (using the abundances of 
ISutherland & Dopitall9 93). Z is the metallicity in solar units, and 
/ovi is the fraction of oxygen atoms in the fifth ionized state. We 
note that the postshock speed is directly related to the postshock 
temperature via 




where we have assumed that the shock is strong (so that v g = v 3 / 4, 
with v s the speed of the shock) and fim p is the mean molecular 
mass. Heck man et"aH <2002h emphasized that for moderately en- 
riched gas, A oc Z (e.g., Sutherland & Dopita 1993). In this fully- 
cooled case, we therefore expect the column density to be indepen- 
dent of the metallicity and physical density. 

We note that ther e is some te ntat ive observational evidence for 
such a picture. Shull et al. 1 2003) and Tumli nson et aljl2005t) each 
found pairs of OVI absorbers displaced by v s > 200 km s _1 from 
associated HI and/or low-ionization state metal absorbers. They 
suggested that the velocity difference may be characteristic of a 
shock. 

Of course, equation Q is only approximate because both t coo \ 
and /ovi depend on temperature. To properly compute the column 
density, we need to trace the (non-equilibr ium) ionic abu ndances 
throughout the cooling (e.g.. iDopita & Sutherland Il996l) . Fortu- 
nately, we can estimate the result reasonably well because /ovi 
peaks sharply around a characteristic temperature, as we show in 
Figure Qj. We have computed the ionization fractions fj on of the 
highly -ionized states of oxygen using Cloudy (version 94. lFerlan3 
2001) for gas in CIE. We see that OVI is extremely sensitive to 
the gas temperature; /ovi ~ 0.3 at T ~ 10 5 ' 45 K but drops 
steeply away from this value. The ionization fraction for OVII, 
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on the other hand, rises rapidly to /ovn ~ 1 in the same tem- 
perature range for which OVI rises, but it remains large for all 
10 5 ' 5 K < T < 10 6 ' 5 K. This is a simple consequence of the 
lone n = 2 electron in OVI, which is much easier to strip than 
the n — 1 electrons in OVII: the ionization potentials of OVI and 
OVII are 138.1 and 739.3 eV, respectively. The plateau in /ovi at 
T ~ 10 6 K and /ovi ~ 0.003 occurs for the same reason. Clearly, 
OVI will come from gas in a narrow temperature range. Figure^ 
shows the ratio of the cooling time to Hq , also computed with 
Cloudy, for a range of metallicities. In the plot we assume that 
A = 10, where A = n/n is the gas density relative to the cos- 
mic mean; of course t coo i oc A" 1 . The cooling time drops rapidly 
as we approach T ~ 10 K from above (partly because of the 
strong OVI AAIOS^ 1038 doublet) and then decreases further as 
OVI recombines. Heckma n et alJ 120021) pointed out that this al- 
lows us to substitute the values of /ovi and t coo \ at the OVI peak 
in equation Q, yielding 

AWi « 10 14 _A cm" 2 (3) 

V 50 km s / 

for enriched gas; note that the result is independent of metal - 
licity and physical density. The column density depends only 

1/2 

on shock velocity (iVovi oc T s ) in this regime. However, 
comparison to the non-equ ilibrium postshock-cooling models of 
DoDita & Sutherland 1 1996) suggests that our prescription under- 
estimates iVovi by about a factor of two for most temperatures. 
This is particularly important for our low-density IGM absorbers: 
the recombination times for the ions and environments of interest 
can be comparable to (or even exceed) the age of the universe. For 
example, at A = 100 and T = 10 6 -10 7 K, we find from Cloudy 
that (for these three ions) the typical radiative recombination times 
are, in most cases, only a few times smaller than the Hubble time. 
We therefore suggest that detailed shock models will be required 
for careful quantitative fits to the observations, especially if some 
shocks occur at lower densities. 

The key assumption of this prescription is that the gas can 
cool pa st T ~ 10 5,45 K . While that is certainly true for most of the 
systems Heckma n et alJ 120021) considered (such as OVI in galactic 
winds), it is only marginally satisfied by IGM absorbers. Even at the 
OVI peak, t coo \ > Hq 1 for Z < 0.05 Zq, and, because metal lines 
dominate the cooling at 10 5 K< T < 10 6 ' 5 K, it increases rapidly 
above T ~ 10 5 ' 5 K as more and more atoms are ionized. Thus, 
tenuous, metal-poor gas and gas shocked to high temperatures will 
not be able to cool through the OVI peak, and the Heckman et al. 
(2002) argument actually only provides an upper limit to column 
densities in the IGM. To account for this possibility, we must limit 
the time over which gas cools to be smaller than Hq . We let 
T m j n = T s [l — (^oicooi) 1 ] be the (approximate) temperature 
to which the gas cools, where the cooling time is evaluated at T s 
and T m i n = if i CO oi < Hq 1 . We then take /ovi to be equal to 
its maximum value in the range (T m i n ,T s ). For equation Q, we 
evaluate the cooling time at this effective temperature. If it exceeds 
the age of the universe, we instead substitute Hq 1 . We note that the 



1 The CIE cooling models of ISutherland & Dopital I 1993ft are often used 
for this purpose. However, those tables modify the [O/Fe] ratio as a func- 
tion of metallicity; because oxygen is particularly important for our pur- 
poses, we have chosen to keep its abundance fixed relative to other metals. 
For consistency with prev ious work, we use the solar-abundance ratios of 
ISutherland & Dopital Il993l) . with the exception of helium (which we set to 
its primordial value for a mass fraction Y = 0.24). 
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Figure 1. (a): Fractional abundance of OVI, OVII, and OVIII (solid, 
dashed, and dotted curves, respectively) if we assume CIE. (b): Ratio of the 
cooling time to the Hubble time, as a function of gas temperature. The solid, 
long-dashed, short-dashed, and dotted curves assume Z = 0.3, 0.1, 0.03, 
and 0.01 Zq, respectively. All curves assume that the gas element has den- 
sity ten times the mean (i.e., A = lOAio, with Aio = 1). 



solution is more uncertain when t coo \ > Hq 1 , because the density 
and metallicity of the gas do affect the total column density. 

We illustrate these features in Figure|2| The thick curves show 
the column density of OVI as a function of shock temperature for a 
range of absorber metallicities and densities. We emphasize that the 
maximum A?ovi is nearly constant when t coo i < Hq 1 . However, 
for small density or metallicity and for large T s , Novi oc AZ. 
Thus if OVI absorbers correspond to shocked regions, we expect a 
characteristic maximum iVovi ~ 10 14 cm" 2 . This is only approx- 
imate for (at least) two reasons. First, the viewing angle relative to 
the shock plane will smear out the sharp maximum. Second, as de- 
scribed above, detailed non-equilibrium calculations show that the 
column densi ties may be so mewhat larger than we predict (by a 
factor of ~ 2: lDopita & Sutherland 1996). The crucial point is that 
the distribution of column densities should depend on the distribu- 
tion of shock temperatures in the IGM. Its maximum is relatively 
robust to the choice of density and metallicity, but the abundance 
of weak absorbers can be quite sensitive to those values. 

Of course, this postshock-cooling model does not just apply 
to OVI. As lHeckman et all 12002) pointed out, it also predicts that 
a range of ionization states should appear and disappear as the gas 
cools. Here we will focus on OVII and OVIII, whose detection 
with X-ray telescopes has received a good deal of attention. Naively 
comparing /ovi and /ovn in Figure^ implies that every OVI 
absorber associated with a shock must also have a large iVovii (al- 
though the converse is not true). We can use the same procedure as 
in equation to estimate the column density as a function of tem- 
perature. The prescription is actually somewhat simpler because 
/ovn is essentially flat. On the other hand, /ovin is large only if 
T a > 10 6 - 5 K, where t coo i > Hq' 1 . Thus we do not expect OVI 
absorbers to host strong associated OVIII absorption, and vice 
versa. Note that this is one important difference between our (time- 
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Figure 2. Column density of cooling material corresponding to a given 
shock temperature. The lower set of thick lines are for OVI, while the up- 
per thin lines are for OVII. Solid, long-dashed, and short-dashed curves 
have A = 10 with Z = 10 -1 , 10" 1 ' 5 , and 10~ 2 Z Q , respectively. Dot- 
dashed and dotted curves have A = 100 with Z = 10" 1 ' 5 and 10~ 2 Zq. 
The top axis shows the halo mass we associate with each infall shock (see 

limited ) model and the fully-cooled predictions of iHeckman et alj 
1 2002). Of course, OVII and OVIII can coexist in the same system. 

The thin curves in Figure [2] show the OVII column densities 
corresponding to each shock temperature. We see that OVII ab- 
sorption should be significant over a much larger range of temper- 
atures than OVI, but because the cooling times are long, the maxi- 
mum Novri is much more sensitive to the gas density and metallic- 
ity. We therefore cannot robustly predict a maximum iVovii. On the 
other hand, note that iVovii is almost always several times larger 
than Novi, even near the OVI peak. 

The postshock-cooling model also predicts physi- 
cal sizes corresponding to each column: t « « 9 t C ooi ~ 
500 (v 9 /50 kms" 1 ) (H tcooi) kpc. (Note that the last factor 
cannot exceed unity in our time-limited picture.) It is possible that 
geometric effects (such as the filamentary structure of the cosmic 
web) could provide more stringent constraints on the sizes (and 
hence column densities) of the absorbers. Also, note that when the 
cooling time is long, the density profile of the gas behind the shock 
can affect the column density (i.e., the assumption of a constant 
density medium may not be appropriate). 



3 THE ABUNDANCE OF OXYGEN ABSORBERS 

If the WHIM gas and the oxygen absorbers are a result of structure 
formation, it is natural to associate them with virialized dark-matter 
halos. We will suppose that each halo has a cross-section to host an 
absorber that is proportion al to the square of its vi rial radius r v i r , 
which we will define as in lBarkana & Loebl 1200 ll) . We can moti- 
vate such a scenario by noting that each halo is surrounded by a 
virial shock that heats accreting material to the virial temperature 
Tvi r . Massive galaxy groups and clusters thus have a hot gas reser- 



Figure 3. Number of virialized objects intersected per unit redshift as a 
function of halo mass. The top axis shows the corresponding virial temper- 
ature. The solid and dotted curves are for z = and z = 0.2, respectively. 

voir within r v h-. With this picture, we can compute the differential 
number dN/dz of absorbers per unit redshift via 

dN f , . , 2 dr 

— = dmn(m)nr viT — , (4) 
dz J dz 

where all quanti ties a re in comoving units. Here n(m) is the 
Sheth & Tormen 1 1999) halo mass function. Figure|3|shows the in- 
tegrand (weighted by mass m). The crucial result is that the num- 
ber of systems intersected is small, of order a few per unit redshift. 
Thus, to reproduce the observed number density of OVI and OVII 
absorbers, we must have absorbers in the IGM outside of virial ob- 
jects; however, the virialized systems could still account for strong 
absorbers, so we cannot neglect them entirely. 

Moreover, we do not expect every halo to host an absorber of 
a given type. As the depth of the potential well increases, the virial 
shock becomes stronger and the gas temperature increases. Thus 
higher ionization states trace gas associated with larger objects, and 
we require some prescription to associate dark-matter halos with 
ion columns. We will now consider two simple possibilities: virial 
shocks and infall shocks. 

3.1 Virialized Halos 

Once gas falls through the virial shock, it is heated to T v i r and be- 
gins to cool radiatively. If t coo i is small, it accretes onto the central 
galaxy, passing through a range of ionization states along t he way. 
If t co o i is large, the gas remains in a diffuse halo IWhite & Frenkl 
1991 ). Virial shocks thus provide an obvious source for hot absorb- 
ing gas. 2 The relevant halo masses can be read off from the top 
axis of Figure|3| which shows T v i r as a function of halo mass: we 

2 We note that both simulations and analytic arguments suggest that virial 
shocks do not form in many systems i Birnboim & Dekel 2003 ; Keres et al. 
2004). However, for the halo masses and redshifts considered here, the 
"hot" accretion phase dominates and virial shocks should be omnipresent. 
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Figure 4. Number of virialized absorption systems with column density larger than some threshold per unit redshift for the OVII and OVIII ions (panels 
a and b, respectively) at z = 0. Long-dashed, solid, and short-dashed curves assume Z = 0.03, 0.1, and 0.3 Zq with A coo i = 200. The dotted curves 
assume Z = 0.1 Zq and A coo i = 50. The thick curves show the total dN/dz, while the thin curves (plotted only for the solid and short-dashed cases) 



show the contribution from systems with t coo i 
everywhere. 



< H Q . In panel b, the dotted curve is coincident with the long-dashed curve and the thin solid curve is zero 



see that massive galaxies and groups correspond to the temperature 
range of interest. 

To compute the absorber st atistics, we will follow the simple 
picture of Wh ite & FrentJ ll99ll) and divide the gas into two com- 
ponents. First, if t coo \ < Hq 1 , we assume that it cools fully. We 
evaluate the cooling time at a density A CO oi (usually 200) reflecting 
the typical density within a halo. In this case, we can compute the 
column density of each system using the postshock cooling model 
described in Sj2| each halo mass then corresponds to a single col- 
umn density independent of the impact parameter (at least if we 
neglect orientation effects). Equation J4} gives the number of ab- 
sorbers per unit redshift provided that we include only those halos 
hosting absorbers of the specified column. From Figure^?, it is ob- 
vious that at virial overdensities, cooling will be complete for halos 
with T vir < 10 6 K. This regime will thus be especially important 
for the OVI absorbers. 



In the other regime, where t coo i > Hq , the gas settles into 
a quasi-equilibrium hot halo with T ~ T v i r . We could again use 
the model of j2|to compute the corresponding columns. However, 
the density structure is likely to be important for virialized halos, 
and reasonably well-motivated models for the profile do exist. We 
will therefore use a simplified version of Perna & Loeb 1 1 99a) . in 
which we integrate along a line of sight through the halo to find 
the column density; obviously the column density does depends on 
the impact parameter in this case. We will assume that a fraction 
fh = 0.5 of an object's baryonic mass is contained in a diffuse 
gaseous halo within r v ir- We will also assume that its density fol- 
lows a beta-model distribution, p oc [1 + (r/Yc) 2 ]^ 1 , where r c 
is the core radius. We will take r c = r v i r /4 = cpr v i r , which 
is approximately consistent with the profiles of luminous clusters 
I J ones &Formanl l999). We will assume isothermal gas, so the col- 



umn density Ni on along a line of sight with impact parameter b is 

2A. c n H {z) r c 



/ion Z [O/H] f h 



m m g{cf)) v /i + ( } 2/ r 2 



xtan 




(5) 



where the metallicity is in solar units, A c is the virialization over- 
density in units of the critical density t Bryan & Norm arTl998l) . and 
g(x) = x 2 (l — x tan _1 a; _1 ). In the slow-cooling limit, we have 
-/Vion oc fhZ. The number of systems per unit redshift is then the 
analog of equation J4}, 



cUV(> N ion ) 
dz 



dr 
dz 



dm n(m) nR 2 (m,Ni 



(6) 



where R(m, TVion) is the maximum impact parameter for a halo of 
mass m within which the column exceeds jVkm. 

Thus, to compute the total number of absorbers from virial- 
ized halos, we begin by dividing halos into fast-cooling and slow- 
cooling regimes (this amounts to finding the halo mass above which 
tcooi exceeds the age of the universe). For masses above this thresh- 
old, we approximate the halo as isothermal and use the halo mass 
profile of equation ^5} to convert to column density as a function of 
impact parameter; the abundance is then given by equation ^6}. We 
need not specify a maximum halo mass, because n(m) falls expo- 
nentially fast beyond the nonlinear mass scale. For smaller masses, 
the column density is determined by the shock physics of equation 
Q independent of impact parameter, and the corresponding cross- 
section is given by equation if we integrate only over those ha- 
los with sufficiently large columns. In this case, the integral con- 
verges because small halos correspond to small virial temperatures 
for which highly-ionized oxygen is vanishingly rare (see Fig.|2j- 

We show our results for OVII and OVIII in Figure [4] The 
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Figure 5. Same as Fig.|4] except for OVI. 



thick curves show the total dN/dz, while the thin curves show 
the contribution from systems with t coo \ < Hq . Long-dashed, 
solid, and short-dashed curves assume Z — 0.03, 0.1, and 0.3 Zq 
with A coo i = 200. The dotted curves assume Z — 0.1 Z© and 
A coo i = 50. As expected from Figure [5] the number of virial- 
ized absorbers is a few per unit redshift, well below the total fre- 
quency found in o bservations iNicastro et all2005h and simulations 
iFane et aljEobl) but consistent with other analytic treatments of 
virialized objects erna & Loeb 1998). Systems with t coo i < H 
can make a significant contribution to the OVII absorbers but are 
unimportant for OVIII . However, they are always relatively weak, 
with a characteristic iVovn ~ 10 15 ' 5 cm' 2 (that is relatively in- 
dependent of metallicity). In contrast, hot halos can contain quite 
large columns, because the path lengths can be long. The abun- 
dance of OVII absorbers is larger than that of OVIII because halos 
near the nonlinear mass scale have T vil < 10 6 5 K, below the ion- 
ization threshold of OVIII. 

Note that the number of high-column-density absorbers in- 
creases if the metallicity decreases. As the metallicity declines, the 
cooling time increases and less-massive halos retain their hot gas. 
This makes up for the decreased column in any one absorber be- 
cause n(m) falls rapidly with mass. This is in many respects an 
artifact of our simplified treatment: for t coo i ~ H^ 1 , we should 
include cooling to compute the detailed column densities. 

Here and below, we present results for z = 0. The statistics 
evolve only modestly to z ~ 1 (e.g., compare the two curves in Fig- 
ure [3). Beyond that, the abundance falls rapidly because the non- 
linear mass scale falls below T vir ~ 10 6 K. 

Figure [5] shows analogous results for OVI. In this case, the 
rapidly-cooling halos are even more important, yielding a charac- 
teristic column density iVovi ~ 10 14 ' 5 cm -2 (see Figure[2j. The 
high-iVovi tail is composed of halos just above the cooling thresh- 
old which still have relatively large /ovi- Thus it only appears if 
the cooling time is long (i.e., Z or A coo i is small) and is subject to 
uncertainties in our simplified cooling model. 

As required by Figure [3] virial shocks underpredict 
d./Vovi/dz and diVovn/dz by factors of several. Clumping will 



further decrease the cross section. If we place the shocked gas into 
rib blobs of radius r < r vlI , only a fraction / cov ~ nb(r /r v i r ) 2 < 
1 of lines of sight through each halo would host absorbers. It could, 
however, increase the column densities. In the slow-cooling case, 
we would have iVj on cx /cov » increasing the amplitude of the 
strong-absorber tail. On the other hand, fast-cooling gas would be 
unaffected. 

Finally, we note that we have allowed the gaseous halos 
in groups to remain relatively dense and hot. It is well-kn own 
that s uch models ov erpredict the soft X-ray background l Pen 
ll99St IWuetaljEoOlh . One solution is to introduce a feedback 
mechanism that decreases t he density of the intragroup medium 
JWu et aljEoOlllXue & Wull2003l) . This would, of course, modify 
our column-density distribution. By slowing the cooling, the prin- 
cipal effect would be to increase the amount of absorption. It could 
also slightly increase the cross section if the gas is spread over 
larger distances. However, simple geometry demands that such pre- 
heating scenarios will certainly not bring the abundance into agree- 
ment with observations; it is more likely to change the column den- 
sities by a modest amount. In principle, the distribution of strong 
IGM absorbers could then be used as a probe of feedback mecha- 
nisms, though the many other uncertainties (such as metallicity and 
non-equilibrium effects) make this a difficult endeavor. 

3.2 Mailing Gas 

To reproduce the large number of absorbers found in simula- 
tions and observations, we must appeal to a population beyond 
virial shocks. This is of course not surprising, because simula- 
tions contain large networks of IGM shocks. To connect these 
shocks with collapsed objects, we will postulate that the temper- 
ature T s h of IGM shocks surrounding each dark-matter halo de- 
pends (exclusively) on the halo mass. One way to motivate such a 
pictur e is if the shocks oc cur during infall onto existing structures. 
Furlanetto & Loeb (2004) used this ansatz, together with spheri- 
cal symmetry, to estimate the distribution of IGM shocks. They ar- 
gued that shocks could plausibly occur when an object breaks off 
from the cosmological expansion (at "turnaround"), because that is 
the point at which the flows begin to converge. Assuming strong 
shocks, the peculiar velocity at turnaround implies T sn ~ T v i r /4. 
The model roughly matches simulation estimates of the mass frac- 
tion and characteristic temperature of the WHIM. Although clearly 
a simplification, this picture does describe many of the major fea- 
tures of the WHIM and provides a well-defined way to convert 
n(m) to the distribution of T sn . Furthermore, spherical infall oc- 
curs on quite large scales and could provide the cross section we 
need. Turnaroun d occurs at r t a ~ 3.3r v h- for self-similar spher- 
ical infall (Bertschinge3 ll98^) in an Einstein-de Sitter universe, 
so the shock networks could be quite large relative to collapsed 
objects. Simulations also show networks of shocks surrounding 
bound structures at a few times the virial radii iKeshet et all2003t 
iNaeai & Kravtsovll2003l) . We will, however, leave the cross sec- 
tion of each infall shock free, because it must depend on the asym- 
metric, filamentary nature of the infall. We will parameterize it as 
X7rr 2 ir . For reference, the top axis of Figure [2] shows T sn (m) in 
this infall model. 

Because of the lack of data pertaining to IGM shocks, a de- 
tailed comparison of this model with observations is difficult. How- 
ever, iRines et all 120031) have used galaxy redshifts to trace the 
infall regions of eight nearby clusters spanning about an order 
of magnitude in mass. Although the galaxies obviously will not 
be subject to shocks, their structure should reflect the underlying 
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Figure 6. Number of absorption systems from infall shocks for the OVI ion 
at z = 0, assuming \ = 8. The long-dashed, thin solid, and short-dashed 
curves assume Z = 0.03, 0.1, and 0.3 Zq, respectively, with A = 10. 
The dot-dashed curve assumes Z = 0.1 Zq with A = 100. The dotted 
curve is the same as the solid curve except that we take T a ^ = T v i r /6. The 
two sets of points show observed results (the solid squares are described 
in the text; the open hexagons are from Danforth & Shull 2005). The thick 
solid curve adds the contribution from virialized systems. 

gas di stribution at least qualitatively. Most importantlv. lRines et alJ 
(2003) found that r- ta ~ 3.5r v i r in all of their clusters. The 
spherically-averaged density profil es fall as ~ r~ 3 -r~ 4 at large 
radii, consistent with simulations I H ernquislll990l : iNavarro et alJ 
Il997l) . This implies that the average density falls by about a factor 
of fifty between the virial radius and turnaround radius. However, 
as we will show in j4]below, the density field is highly inhomoge- 
neous within this region. For simplicit y, we will assume a constant 
density behind these shocks. Finally, iRines et all 120031) compute 
the velocity dispersion profiles. Between r v i r and rt a , the disper- 
sion typically declines by a factor of two or so. This is nicely con- 
sistent with our assumption that T s h ~ TU r /4. 

Given this prescription, we can compute the expected number 
of absorbers using the postshock cooling model of ij2| Here we will 
use this model for both fast-cooling and slow-cooling shocks, be- 
cause (unlike for virialized halos) we have no well-motivated model 
for the density profile of IGM filaments. However, we caution the 
reader that geometry (such as the thickness of filaments) could still 
ultimately determine the column densities. We also note that, for 
these modest densities A ~ 10, the recombination times can ap- 
proach or even exceed the age of the universe, so our assumption 
of local CIE in calculating the abundances may be problematic (see 
discussion in Sj2j. To compute dn/dz, we use equation J4j, with 
the integration range set by the column density corresponding to 
each halo mass. Operationally, the calculation is similar to that de- 
scribed in EI3.1I except that equation Q applies to both fast and 
slow-cooling shocks. 

Given the simplifications of our model, a detailed compari- 
son to observations is unwarranted. However, we do wish to esti- 
mate the required \ as we H as check whether the derived column 
densities are reasonable. As such, we begin with OVI, which is 



by far the best studied ion. The filled squares in Figure |6| show 
a dis tribution compiled from several lines of sight in the litera- 
ture iTripD et alJbOOCt fSembach et ail 1200 it iRkhter et aljEbol 
IProchaska et alj|2004) : these lines of sight have the best-defined 
samples. 3 We have included all claimed detections (even those 
labeled "tentative" by each survey's authors). The error bars as- 
sume Poisson statisti cs. The open hexago ns show the more sys- 
tematic survey of Danforth & Shull 1 2005), who searched for OVI 
absorbers in 129 known Lya absorbers along 31 lines of sight stud- 
ied with the FUSE satellite. We refer the reader to that paper for 
a detailed discussion of the errors. Note that the t wo data sets are 
reasonably consistent, though the Danforth & Shull (2005) data do 
lie below our compilation by a small amount. 

The curves in Figure|6|show the results from our model assum- 
ing x = 8. The long-dashed, solid, and short-dashed curves assume 
Z = 0.03, 0.1, and 0.3 Z Q , respectively, with A = 10. The dot- 
dashed curve assumes Z = 0.1 Zq with A = 100. The dotted 
curve is the same as the solid curve except that we take somewhat 
weaker shocks, with T s h = T v i r /6. The thick solid curve adds the 
contribution from virialized systems for a better comparison with 
the observations. 

In all cases, the number of OVI absorbers falls sharply at some 
characteristic column density. This is simply the maximum allowed 
for fast-cooling shocks. In contrast to virial shocks, however, the 
densities are typically small enough that the maximum column den- 
sity does depend on the metallicity and physical density (unless 
Z ~ 0.3 Zq or A ~ 100, in which case the cooling times be- 
gin to fall below the age of the universe). Below the characteristic 
maximum, the number of absorbers rises rapidly because /ovi is 
so sensitive to the shock temperature. We thus expect infall shocks 
to contribute only to the range TVovi ^ 10 14 cm ~ 2 - Note that the 
sharp rise in the dot-dashed curve at JVovi ~ 2 x 10 13 cm -2 oc- 
curs because of the plateau in /ovi at T ~ 10 6 K (see Fig.0. Of 
course, the maximum Vovi does not decrease for T s h = T v i r /6 
because that has no effect on the cooling. It decreases dN/dz by a 
small amount because it shifts the halo mass corresponding to each 

Tsh- 

Thus, if all the known absorbers are to be attributed to shocks 
(either infall or virial), we find that \ ~ 8; i.e., the infall shocks 
subtend an area eight times larger than the halo itself. While large, 
this is smaller than the expectation from pure spherical infall so 
may not be unreasonable. Moreover, as we will discuss in ^be- 
low, a number of the observed systems appear to be photoionized 
and have a different origin, which would allow the covering factor 
of infalling gas to decrease by about a factor of two. This would 
decrease the thin curves in Figure |6| by the same factor: note that 
we would still expect a sharp decline near the characteristic maxi- 
mum column density of the shocks. Because stronger systems only 
arise in virial shocks, the covering fraction within virialized objects 
can be separately calibrated from the abundance of strong absorbers 
(though of course that could also be a function of halo mass). In the 
simplest model, a significant population of photoionized absorbers 
would fix x ~ 4, allowing IGM absorbers to be more filamentary, 
while leaving the virial shock spherical. In any case, we find it re- 
assuring that a reasonable choice for \ can simultaneously match 
the observed abundance and column-density distribution, at least 
roughly. 



3 However, we have not attempted to account for the different detection 
limits of each survey. As such, the statistics at ATovi ~ 10 13 cm -2 may 
not be reliable. 
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Figure 7. Same as Fig.|^] except for the OVII and OVIII ions (panels a and b, respectively). Note that the ordinates have different scales. 



We repeat the same exercise for OVII and OVIII in Fig- 
ure Q We have again fixed x = 8. First, note that the number 
and characteristic columns of OVII absorbers are much larger than 
for OVI absorbers. This follows easily from Figure [2] because 
OVII dominates over a much larger temperature range (includ- 
ing the OVI peak). On the other hand, the maximum iVovn is 
poorly constrained: because the cooling times always exceed the 
age of the universe, the cutoff is proportional to metallicity and 
density. Interestingly, OVII drops off steeply at this maximum 
and is nearly flat below it, so most absorbers of a given metal- 
licity and density should cluster around one particular column. 
This occurs because every system cools over H^ 1 and /ovn is 
nearly independent of temperature. In reality, the cutoff will be 
smeared by the range of metallicity and density as well as ge- 
ometry. For X-ray lines, the comparison to observations is much 
more difficult given current technology. There is only one well- 
studied line of sight (with secure detections of two absorbers); for- 
mally this implies d N/dz(N vu > 7 x 10 14 cm" 2 ) = 67±|| 
iNicastro et alj|20o3) . Our model is well within this range. More 
interesting is to note that the observed columns (7 x 10 14 and 
10 15 cm -2 ) are both near the cutoff and appear to require A > 30 
or Z ~ 0.3 Zg, still well within reasonable bounds. Simulations 
suggest dA7dz(i Vovii > 7 x 10 14 cm" 2 ) ~ 10-30 , close to out- 
model with x = 8 <Fang et alj200llChen et all2003h . They do not 
show as severe a break as our model, probably because of scatter in 
the metallicity and geometric effects. 



FigureQ? shows that the expected number of OVIII absorbers 
is about half that of OVII, although the column densities can be 
comparable. OVIII is rarer simply because the characteristic shock 
temperature is < 10 6 ' 5 K iFurlanetto & Loebl2004) . so there is not 
as much gas with large /oviii- As a result, virialized systems make 
a larger relative contribution. Also, note that OVIII does not show a 
sharp cutoff; this is because /oviii varies rapidly with temperature. 



4 COMPARISON TO SIMULATIONS 

We have shown that the spherical-infall model described in the pre- 
vious section can explain many of the key properties of the ab- 
sorbers. However, it also makes a number of important assumptions 
and simplifications. Here we compare our model to simulations in 
order to highlight its strengths and weaknesses. 



4.1 The Cosmological Simulation 

We use the z — dark-matter, stellar-particle, gas-temperature 
and density outputs of a large scale L = 25fo -1 Mpc cosmologi- 
cal hydrodynamical simulation, with 768 3 fluid elements and 384 3 
dark-matter particles (each with mass 2x_10 7 hr^Ma). The code 
structure is similar to that in|C^n_&_Ostrikej ll992albl) with some 
significant changes l Cen & Ostriker 2000). It is described in detail 
in lNaeamine et alj i200l|). The c hosen cosmology is the "concor- 
dance model' ' (Wang etalJl2000h . a flat low-density (fl m = 0.3) 
universe with a cosmological constant (S1a = 0.7), and h — 0.67. 
The baryon density (originally £lt = 0.049) has been scaled to 
fit = 0.046 to match the value used in our analytical model. The 
galaxy pa rticles a re formed through the recipe described in the Ap- 
pendix of iNaeamine et all fcOOOf) . Each stellar particle is assigned 
a position, velocity, formation time, mass, and metallicity at birth. 
The stellar particle is placed at the center of the cell with a velocity 
equal to the mean velocity of the gas in the cell. It is then followed 
by a particle-mesh code as a collisionless particle interacting grav- 
itationally with the dark matter and the gas. We use the galaxy and 
group position information obtained from these stellar particles by 
INaeamine et alJfcOOll) . 

We will primarily be concerned with the IGM around bound 
objects. We use the dark-matter halo pa rameters for this simula- 
tion obtained by Nagamine et alj j2 00ll) with the HOP grouping 
algorithm of Eisenstein & Hu 1 1998). We then search for the ha- 
los that are most closely coincident (AL < 68 h~ x kpc) with the 
galaxy-particle groups and use the properties of these dark-matter 
halos. The virial radius i s obta ined from halo mass using equation 
(24) of Bar kana & LoeblfeoOll) . Around the center of each halo, we 
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obtain the volume-averaged radial distributions of matter, temper- 
ature, and metals in the inner regions abutting the galaxy-particle 
locations in the z = output: 

Jd 3 xp(x)Q(x) 
J d 3 xp(x) 

where Q(x) is the parameter value (A, T, /j OJ1 , etc.) in the volume 
element d 3 x. For the sake of transparency, we do not directly use 
the metallicities reported by the simulation. Instead we consider 
two metallicity prescriptions. The first scales metallicity with den- 
sity as described bv lCroft et all 1200 lh . Z = 0.005 Z x yjn/n, 
with a maximum metallicity of Z — 0.33 Zq. The second pre- 
scription simply uses a constant metallicity Z = 0.1 Zq. In or- 
der to study the OVI , OVII , and OVIII abundance distributions, 
we compute the ionization fraction as a function of temperature 
as illustrated in Figure With f\ orL for each cell, plus the den- 
sity and metallicity, we compute the total ionic abundance relative 
to hydrogen. We scale the distributions by the virial radius corre- 
sponding to each halo and follow them out to five virial radii. We 
obtain the volume-averaged distributions for several virial-radius 
bins, r vir = 240-300, 360-480, and 600-720 kpc; these corre- 
spond to roughly (1-2), (3-7), and (14-24) x 10 12 M©. 

4.2 Absorbers Around Virialized Halos 

The key ansatz in our infall-shock model is that the characteris- 
tic IGM temperature surrounding a collapsed object increases with 
its mass, even on scales much larger than r v i r . Figure [8] shows 
spherically-averaged temperature profiles around halos in the sim- 
ulation (normalized by r v i r ). 4 We see that isolated galaxies tend 
to lie in cooler regions of the IGM than massive groups, even at 
several times r v i r . The trend appears to break down at r > 3r v j r 
for the smallest objects. This is in a temperature range where cool- 
ing can be significant, which may affect the results. Also note that 
there is no characteristic radius at which the temperature increases, 
as would be expected for pure spherical infall. Thus the shocks ap- 
pear to be complex and distributed, with a range of temperatures. 
The typical temperature is, however, a few times smaller than T v i r , 
consistent with our model. 

Figure|5|shows that it may be more accurate to associate dif- 
ferent oxygen ions with different halo masses. The three panels 
show the spherically-averaged ionization fractions for each of our 
ions as a function of normalized radius. The curves correspond to 
the same mass ranges shown in Figure|8| We see that OVI tends to 
be associated with the smallest-mass objects, although even these 
halos are too hot to sustain a large /ovi inside r v i r . On the other 
hand, OVIII is most prevalent around the highest-mass objects (al- 
though not within the virialized regions, again because the gas is 
too hot) and is extremely rare around the low-mass objects. Finally, 
OVII appears in all systems because it dominates over such a wide 
temperature regime; only the virialized gas in extremely massive 
objects is hot enough to destroy OVII . 

Thus, we can be confident in expecting highly-ionized ab- 
sorbers to lie near massive collapsed objects. However, the simula- 
tions also show that a simplified spherical model is not accurate. In 
Figure fTol we have computed dn/dz for each ionic absorber using 
our spherical profiles. We have used equation |6j with the column 
densities determined by integrating through the average profiles. 

4 At r -C r v ; r , the temperature is strongly affected by star formation and 
feedback, so we will ignore that portion. 
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Figure 8. Spherically-averaged temperature profiles around virialized halos 
in a cosmological simulation, where all distances are scaled to the virial 
radius. The three curves show different mass ranges. 
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Figure 9. Spherically-averaged ionization fractions around virialized halos 
in a cosmological simulation, where all distances are scaled to the virial ra- 
dius. The top, middle, and bottom panels are for the OVI, OVII, and OVIII 
ions. In each panel, the three curves show different mass ranges, as in Fig.l8l 

In each case we truncate the profile at 5r v ir. The dotted, dashed, 
and solid curves show the expected statistics for OVI, OVII, and 
OVIII, respectively. First consider the thin curves, which assume 
a constant Z = 0.1 Z© metallicity (as in our fiducial analytic 
model). We see dn/dz ~ 65 for each ion at sufficiently small 
column densities. This is an artifact of truncating the integrating at 
5r v i r . The profiles in Figure|5|tend to approach nearly constant val- 
ues at large radii, so for large impact parameters the column density 
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Figure 10. Absorber statistics for spherically-averaged profiles taken from 
the simulations. The different line styles show the three different ions. The 
thin curves assume a uniform Z = 0.1 Zq metallicity; the thick curves use 
a density-dependent metallicity prescription (see text). 



is essentially proportional to the maximum radius; only well above 
this turnover can the column densities be trusted. 

In that range, OVII and OVIII do not present dramatically 
different pictures from the analytic model. Because /ovu is nearly 
temperature-independent over a broad range, we might expect 
spherical averages to do a passable job with that ion. Indeed, the 
number density does agree roug hly with our m odel and with more 
detailed simulation analyses iFang et al]l2002t IChen et alJl2003h . 
But it is certainly not so good that we can claim support for our 
model. OVIII is considerably worse, and the OVI results do not 
appear to have converged and disagree strongly with our model. 
This is also not surprising, because OVI is often in a fast-cooling 
regime and is so temperature sensitive that averaging over large 
volumes artificially suppresses its importance. Ray tracing appears 
vital to predict the OVI absorber abundance, as would be expected 
in a shock-cooling model. 

The thick curves in Figure fTOl assume the density-dependent 
metallicity distribution described above. The number of absorbers 
declines sharply in this scenario: at large radii the mean physi- 
cal densities are near the cosmic mean, so the overall metallic- 
ity is small. As a result, high column densities require small im- 
pact parameters. One fortunate consequence is that the curves con- 
verge better. However, much more significantly, the number of ab- 
sorbers isf^below the resul t of more complete simulation analyses 
<Cen et all200ltlFang & Brvanl200ltlFang et alJioollChen etai] 
2003). Thus, the spherical averages do not accurately describe the 
absorbers in simulations. Of course, the reason is the dumpiness of 
the IGM. 

Figure fTTI shows phase diagrams of the IGM gas surround- 
ing halos with M — 3-7 x 1O 12 M0 in four radial slices. The 
np yi, npyn, and no ym-averaged contours were obtained using 
the ICroft et all 1200 ll) prescription for the metallicity. Hotter gas 
tends to sit at higher densities (and hence occupy a small volume). 
This correlation holds even at A ~ 1-10, well below virial over- 



densities. Moreover, the mean temperature and density decrease 
and the scatter in the relation increases as one moves away from 
the central galaxy. However, regardless of radius there is a non- 
negligible fraction of gas at high temperatures comparable to T v i r . 
This gas (which also has a relatively large density) is presum- 
ably responsible for most of the oxygen absorption. Thus, as can 
be seen from Figure ITTb . the higher-ionization states will tend to 
probe the densest filaments (and the most enriched gas, on aver- 
age) around each structure, in addition to the largest structures. 
Meanwhile most of the volume is filled with cooler, low-density 
gas that has a small metallicity given our prescription. The averag- 
ing procedure smooths the high-density gas, boosting the column 
density along all lines of sight under a constant metallicity. But 
when the metallicity is density-dependent, the tenuous gas is so 
metal-poor that the dense clumps essentially disappear during av- 
eraging. These clumps are reflected in the fragmentary nature of 
the density-averaged contours at higher A. As one moves outward, 
they represent an ever decreasing fraction of the volume but con- 
tinue to contribute significantly to novll, roovin, and even novi- 
The distribution is highly complex: there is a wide range of temper- 
atures for any given density (spanning about a decade for A ~ 10), 
which reflects the non-uniformity of the infall process. 

How can we reconcile the apparent success of our analytic 
treatment with the clear failure of spherical averaging in simula- 
tions? The key reason is that, with the exception of hot virialized 
halos, we did not use spherical geometry in the calculation; instead 
we used only the infall velocities and a cooling time (capped at 
Hq ). The physics of shock cooling then predicted the column 
densities without recourse to the gas distribution. Thus we only 
need modify our picture to describe a network of inhomogeneous 
shocks with total cross section nxr^ iT rather than a single coherent 
structure. Of course, it is still an open question whether our model 
provides an accurate description of these networks even on a mean 
level (specifically the shock temperature). In the future, quantitative 
treatments must also consider the importance of the scatter seen in 
Figure fTTI 



5 PHOTOIONIZED GAS 

To this point, we have focused on systems that are cooling after 
being collisionally ionized by a shock. In simulations, such sys- 
tems make up only about half of the OVI absorbers ICen et aU 
l200ltlFang & Brvan 2001). The others are photoionized; such sys- 
tems are rare for the strong absorbers but dominate at iVovi ~ 
10 13 cm~ 2 . In the simulations, these absorbers are typically cooler 
and lower density than collisionally-ionized systems. Modeling of 
the observed O VI systems also suggests that about half are pho- 
toionized (e.g., IProchaska et alj |2004). However, this interpreta- 
tion is subject to uncertainties in the relative metal abundances, 
the ionizing background, and whether the medium is multiphase. 
Reliably separating this population is crucial because it typically 
has T ~ 10 4 K, which means that these systems are not part of 
the canonical WHIM. We have shown in <)2|(see Fig.|2j that post- 
shock cooling predicts that individual OVI absorbers should have 
OVII columns about an order of magnitude larger than iVovi- This 
is larger than /ovi//ovii in CIE at the OVI peak because most 
shocks begin at higher temperatures and cool into the OVI range, 
so OVII probes a longer path length through the postshock region. 
We now calculate the associated OVII columns if photoionization 

dominates. 

We have again used Cloudy (version 94, iFerlandboOll) to 
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Figure 11. Contour plots, in the density-temperature plane, of the volume-weighted (grey-scale), novi _we ighted (dashed), novn-weighted (solid), and 
"OVlll-wdghted (dotted) distribution of gas within 0-lr v ; r (panel a), l-3r v j r (panel b), 3-5r v ; r (panel c), and 0-5r v j r (panel d) of halos with r v ; r = 360- 
480 kpc (or M = 3-7 X 1O 12 M ). The contours enclose 0.95, 0.75, 0.5 and 0.25 of the total volume or ion number. 



compute /ovi and /ovn assuming ionization equilibrium with an 
input radiation background. We constructed a grid in density and 
temperature with spacings A log n^r = 0.5 for —7 < log n/j < 
-3 (in units of cm" 3 ) a nd AlogT 1 = 0.1 for 3 < logT < 7 
(in K). We use the z = O lHaardt & Madaul feOOll) ionizing back- 
ground normalized to a total ionizing rate T_i2 = 0.084 (in units 
of 10~ 12 s -1 ) . This roughly matches the observational estimate of 
iDave & TripclfcOOll) . 

Figure fT2l summarizes our results. Panel (a) shows the log- 
arithm of the ratio /ovi//ovii- Dashed (solid) contours indicate 
/ovi//ovii > 1 (< 1). To delineate those regions of phase space 
where we could expect OVI absorbers in the IGM, panel (b) shows 
the logarithm of the path length (in kpc) that would be required to 
produce an absorber with iVovi = 10 13 cm -2 if Z — 0.1 Zq. 
Systems with L < 100 kpc are confined to gas with T < 10 5 5 K 
and tih ~ 10~ 6 -10 -4 cm -3 or to relatively dense gas lying 
close to the ionization temperature of OVI. In the former case, 
we find that /0V1//0V1I > 1, unless the gas has T ~ 10 5 ' 45 K 
(in which case collisional processes dominate and it is part of the 
WHIM) or is extremely tenuous and spatially large. Thus, pho- 
toionized absorbers have smaller associated X-ray absorption than 
in the postshock cooling model and comparison of /ovi//ovii of- 
fers a powerful test of the nature of the absorbing gas. Upper lim- 
its of -/Vovn ~ 3iVovi can robustly identify those absorbers with 
T < 10 5 K and separate WHIM systems from cooler photoionized 
gas. 

This conclusion holds even if the ionizing background has a 
different amplitude or shape. The ratio /ovi//ovn is nearly in- 
dependent of r_i2. The required path length for low-density cool 



gas increases slightly as T_i2 increases, making tenuous absorbers 
more unlikely. The shape of the metagalactic background matters 
little since the ionization potentials of OV and OVI are so closely 
spaced (1 14.2 eV and 138.1 eV, respectively). We would have to in- 
troduce a (physically unmotivated) sharp feature to the background 
in this energy range to substantially change the ion ratios. We have 
verified this with the Haardt & Madau 1 1996) ionizing background, 
which includes only q uasars and is significantly harder than the 
lHaardt & Madaul 1200 ll) background (which also includes galax- 
ies). The corresponding Figure fT2l looks identical. 

It is worth emphasizing that /ovi//ovii is less subject to 
modeling uncertainties than constraints that use only UV lines (e.g., 
IProchaska et alJ l2004) because it depends only on a single atomic 
species. We also note that other ionization states could be use- 
ful. For example, shocks that produce OVIII have t coo i 3> Hq 
and should not contain substantial OVI . Thus a large measured 
■fo vi//oviii would help to identify photoionized systems. Note 
that Heck man et alJ j2002) imply that large columns of OVI and 
OVIII can coexist in the same absorber. This is because they im- 
plicitly assume that the material cools fully; while reasonable for 
Galactic syst ems, this limit is not relevant for IGM absorbers. 

Finally, iHeckman et alJ J2002t) also pointed out that the 
linewidth — column-density relation offers another way to discrim- 
inate between photoionization and postshock cooling. 

6 DISCUSSION 

In this paper, we have attempted to elucidate the physical ori- 
gin and properties of IGM oxygen absorbers (OVI, OVII, and 
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Figure 12. (a): Contours show /ovi//ovil as a function of gas den- 
sity and tempera t ure, a ssuming photoionization equilibrium in a z = 
lHaardt & Madaul i200lh ionizing background with T_i 2 = 0.084. The 
dashed contours indicate /ovi//ovil > 1 while the solid indicate 
/ovi//ovii < 1; they are labeled with the logarithm of the ratio, (b): 
The path length required for a total column density TVovi = 10 13 cm -2 if 
the absorber has Z = 0.1 Zq. The contours are labeled with the logarithm 
of the distance in kpc. 



OVIII) with a set of simple analytic arguments. While simulations 
appear to predict approximately the correct number of abs orbers 
ICen et alfeOOltlFang & Brvarl200lUFang et all2002l:IChen et alJ 
120031) . such studies have only examined the physical origin of the 
absorbing gas in a cursory fashion. We first argued that the typical 
maximum column density of OVI absorbers (Nqvi ~ 10 14 cm~ 2 ) 
follows naturally from postshock cooling lHeck manet"afl l2002). 
suggesting that at least some absorbers are generated by structure- 
formation shocks. We then argued that, if these shocks are associ- 
ated with bound structures, we must have a (more numerous) pop- 
ulation of significantly hotter shocks around larger objects. These 
shocks can be traced by OVII and OVIII; the former is more 
useful because it dominates over such a broad temperature range. 
Moreover, shocked OVI systems must all have strong OVII ab- 
sorption. This contrasts with photoionized OVI absorbers, which 
should have OVII columns not much larger than the OVI columns. 

We showed that a simple picture with virial shocks and "in- 
fair' shocks can account for the observed number density of ab- 
sorbers if each bound halo has a total cross section several times 
larger than 7r?- 2 ir . The inferred cross section depends on the frac- 
tion of observed OVI systems that are photoionized; if it is neg- 
ligible, we find x ~ 8. If about half of the weak absor bers are 
photo ionized (as suggested by some recent data; Prochaska et al. 
2004), x ~ 4; however, this does not necessarily affect the pre- 
dicted abundance of strong absorbers, which mostly arise from 
virial shocks. In our picture, the temp erature of the infall shoc k 
increases with the central object's mass iFurlanetto & Loebl2 004 ) . 
Comparison to numerical simulations shows that these shocks are 
complex, with dense clumps of metal-enriched shocked gas embed- 
ded in a tenuous weakly-absorbing medium. They do confirm our 
key assumption, the association of hotter shocks with more massive 
systems (albeit with substantial scatter). In the future, we hope that 
more detailed comparisons to simulations will test other elements 
of our picture. 

One of our most interesting results is that, while OVII and 
OVI can be associated with the same shocked system (even ne- 
glecting multiphase media from thermal instabilities), they do not 



necessarily sample the same physical volume. As the gas flows 
away from the shock, it cools over a range of temperatures. OVII 
is presumably associated with gas closer to the shock, while OVI 
probes more distant gas. This calls into question modeling of 
the absorbers with single-temperature CIE gas. Figure [2] shows 
that gas with a more or less constant iVovii can correspond to a 
wide range of TVovi because of the range of temperatures through 
which the gas element cools, even if it remains near equilibrium 
at a ny single temperature (wh ich may not be a valid assump- 
tion; |Do^ita&_^utherlanfJl996). The crucial point is that gas with 
A ~ 10-100 and T ~ 10 5 -10 7 K has a cooling time comparable 
to (though often a few times larger than) the Hubble time, so the 
temperature varies across the system. 

The association of oxygen absorbers with virialization and in- 
fall shocks obviously implies that these systems must be spatially 
correlated with galaxy groups and clusters. This is not unique to 
our picture; any model that associates absorbers with large-scale 
shocks in the cosmic web must show a similar correlation. There 
appears to be a relatively strong correlation between OVI absorbers 
and galaxy groups: most lines of sight s how gal axy overdensi- 
ties coincident with the absorber locations l TrioD & Savag ell200ri 
Savage et al. 2002; Richteretal. 2004), although it has not been 
statistically quantified and the correspondence is by no means per- 
fect. Typical separations are < 1 Mpc, suggesting that the correla- 
tion must extend significantly beyond the virial radius. Perhaps the 
clearest example lies along the line of sight toward PK S 2155-304 , 
which has a pair of intergalactic OVI absorbers l Shu ll et all2003l) 
offset by ±400 k m s~ from a galaxy group. Using Chandra, 



Fang et al. 1 2002) also found an OVIII absorber coincident with 
the galaxy group, although subsequent observations with XMM 
failed to confirm this detection. Because the line of sight passes 
near the barycenter of the group, the OVIII could correspond to 
virialized gas while the OVI systems probe gas that has not yet 
reached the virial radius but has nevertheless been shocked during 
infall. This is somewhat similar to the interpretation of Shu ll et alJ 
(2003), who argued that the OVI and OVIII could trace multi- 
phase infalling gas, although in their case they placed the gas in 
a single shocked regio n. However, w e note that the OVII absorbers 
discovered bv lNicastro et alJ 120051) appear to show no correlation 
with galaxy systems. A more detailed study of the relation between 
galaxies and the known absorbers should shed light on the physical 
origin and significance of these systems. For example, photoion- 
ized OVI absorbers need not correlate with galaxies because they 
are not part of the filamentary structure. 

Our model may also apply to the Local Group absorber. 
iNicastro et alJ j2002) noted that all three ionization states appear 
along the line of sight to PKS 2155-304 with AT vi ~ 10 14 citT 2 , 



Novn 



4 x 10 15 cm 



and iVovin 



5 x 10 15 cm 



They modeled the absorber as a single phase and favored a model 
with low-density photoionized gas. The comparable abundance of 
OVIII and OVII contrasted with the relatively large OVI column 
ruled out CIE. As a result the absorber was forced to have a physical 
size ~ 5 Mpc. In our picture, the transitions sample different spatial 
locations and shocks. We can easily accommodate Novri/Novi at 
infall shocks associated with halo masses M ~ 10 13 M Q . The col- 
umn density iVoviii would be too small in such a shock; however, 
our privileged position near the center of the Milky Way and Lo- 
cal Group implies that this transition could be sampling a second, 
nearby shocked component at T v i r ~ 2 x 10 6 K. Thus the X-ray 
absorbers could be cosmological but no more than a few hundred 
kpc in extent. This is reassuring because, if every dark-matter halo 
comparable to the Local Group had a 5 Mpc hot gaseous halo, we 
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would expect to find several times as many high column-density 
OVII and OVIII absorbers as we see or as appear in simula- 
tions <Fang et al][2002h . In that case, the PKS 2155-304 line of 
sight would have to be unusual; however, comparable absorption 
colum ns appear a long many other lines of sight through the Local 
Group I Rasmussen et al. 20031). 

Our simple model is clearly not sophisticated enough for de- 
tailed predictions of the absorber statistics. For those purpos es, cos- 
mological simulatio ns are much more appropriate I Cen et al. 2001 ; 
iFang & Brvanl200ltlFang et alj2002h . fchen et aljbo03l) have pre- 
sented the most detailed treatment to date. They showed that mod- 
ern simulations do reproduce the observed number densities, at 
least to the (admittedly limited) accuracy of the current measure- 
ments and given an ad hoc (though reasonable) metallicity distri- 
bution. We have taken a complementary approach using analytic 
arguments to estimate the typical column densities and to param- 
eterize the number densities in terms of the known abundance of 
collapsed objects. We believe that a combination of the two ap- 
proaches should yield a much deeper physical understanding of the 
substantial baryon reservoir probed by these observations. 

We are grateful to R. Cen and J. R Ostriker for the use of 
their simulation results in Section|4]and to C. Danforth and J. M. 
Shull for sharing the data used in Fig. [6] We thank T. Fang and R. 
Dave for helpful conversations, A. Kravtsov for enlightening dis- 
cussions, and M. Sako for clarifying issues associated with some of 
the early reports of IGM X-ray-absorption detections. SRF thanks 
the Kavli Institute for Theoretical Physics, where part of this work 
was completed under NSF PHY99-07949. This work was sup- 
ported at Caltech in part by NASA NAG5-11985 and DoE DE- 
FG03-92-ER40701. 
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